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We address the constraints on the SUSY seesaw parameters arising from Lepton Flavour Violation observ- 
ables. Working in the Constrained Minimal Supersymmetric Standard Model extended by three right-handed 
(s)neutrinos, we study the predictions for the branching ratios of lj — » l^-y and lj — » 3£j channels. We impose 
compatibility with neutrino data, electric dipole moment bounds, and further require a successful baryon asym- 
metry of the Universe (via thermal leptogenesis). We emphasise the interesting interplay between 813 and the 
LFV muon decays, pointing out the hints on the SUSY seesaw parameters that can arise from measurements 
of #13 and LFV branching ratios. This is a brief summary of the work of Ref . Q. 



1. Introduction 

Supersymmetric (SUSY) extensions of the Standard 
Model (SM), including three right-handed neutrino 
superfields, are well motivated models which can ac- 
commodate a seesaw mechanism [2;], and at the same 
time stabilise the hierarchy between the scale of new 
physics and the electroweak (EW) scale. One of the 
most striking phcnomenological implications of SUSY 
seesaw models is the prediction of sizable rates for 
lepton flavour violating (LFV) processes Q , many or- 
ders of magnitude larger than those expected from the 
SM seesaw. In this sense, the lj — > h 7 and lj — > 2>U 
(i 7^ j) lepton decay channels, as well as ijl — e conver- 
sion in heavy nuclei, are among the most interesting 
processes. Experimentally, the most promising decay 
is the fj, — > e 7 process, which exhibits the most strin- 
gent present bounds, and offers a significant improve- 
ment regarding the future sensitivity. 

Given the fact that both light and heavy neutrinos 
enter in the determination of the LFV rates (via the 
Yukawa interactions), a powerful link between the low- 
and high-energy neutrino parameters can be obtained 
from these LFV processes. From the requirement of 
compatibility with current LFV bounds and with low- 
energy neutrino data, one can then extract informa- 
tion on the heavy neutrino sector, thus providing an 
indirect access to the heavy neutrino parameters. 

In Ref. [l| , we have systematically explored the sen- 
sitivity of LFV processes to 6*13 in a broad class of 
SUSY seesaw scenarios, with different possibilities for 
the mixing in the neutrino sector. We have also in- 
corporated in our analysis the requirement of gener- 
ating a successful baryon asymmetry of the Universe 
(BAU) via thermal leptogenesis Q. In particular, we 
have shown that various of the lj —> and lj — > 3 U 
(i j) channels indeed offer interesting expectations 
regarding the sensitivity to #13. This sensitivity to 
6*i3 had been previously pointed out @ , for some 
specific seesaw cases. 



Ultimately, and as shown in [l[ , the impact of a po- 
tential #13 measurement on the LFV branching ratios, 
together with the current and future experimental 
bounds (measurements) on the latter ratios, may lead 
to a better knowledge (determination) of the heavy 
neutrino parameters. 



2. LFV within the SUSY Seesaw 

The leptonic superpotential containing the relevant 
terms to describe a type-I SUSY seesaw is given by 
W = N c Y v LH 2 + E c YiLH l + \N c m N N c , where 

N c is the additional superfield that contains the right- 
handed neutrinos and their scalar partners, Y\^ v are 
the lepton Yukawa couplings and to at is Majorana 
mass. Henceforth we will assume that we are in a 
basis where Yi and Tom are diagonal in flavour space. 
After EW symmetry breaking, the full 6x6 neutrino 
mass matrix is given in terms of the 3x3 Majorana 
mass matrix tojv 5 and the 3x3 Dirac mass matrix 
itid = Y V V2, where Y v denotes the neutrino Yukawa 
couplings and 1^1(2) are the vacuum expectation values 
of the neutral Higgs scalars, with v\[2) — v cos(sin)/3 
(v = 174 GeV). 

In the seesaw limit, v <C tojv, we obtain the see- 
saw equation for the light neutrino masses, m v = 
—m^m'^mD ■ The diagonalisation of the full neu- 
trino mass matrix leads to the six physical Ma- 
jorana states: three light vi and three heavy 



states Ni. Their masses are given by to 

,) and m 



C/ MNS m ^ C/ MNS = diag(?7i ;/ 



diag 
/ 

diag 



diag (to^ , toat 2 , toat 3 ). We use the standard parame- 
terisation for the Maki-Nakagawa-Sakata unitary ma- 
trix £/mns |Z|j m terms of three mixing angles 812, 6*23 
and 6*i3, and three CP violating phases, 5, <j)\ and (j>2- 
Following the parameterisation proposed in 0] , the 
solution to the seesaw equation can be written as 



m D 



m d ™ s R 



^MNS ' 



(1) 
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where R is a generic complex orthogonal 3x3 ma- 
trix, defined by three complex angles Oi. This param- 
eterisation allows to accommodate the experimental 
data, while leaving room for extra neutrino mixings, 
in addition to those in C/mns- It further shows how 
large Yukawa couplings Y v ~ 0(1) can be obtained by 
choosing large entries in m^ ag . 

In our analysis, we have considered scenarios of hi- 
erarchical heavy and light neutrinos, mjv x <S mN 2 ^ 
m,N 3 and m Ul <C m v . 2 <ti m U3 , with m^ 2 
andm^ = Am 2 atm +r 



Regarding the numerical es- 



timates, we have used Ato 2 o1 = 8x10 5 eV 2 , Am 2 tm = 
2.5 x 10~ 3 eV 2 , 6»i2 = 30°, 23 = 45°, 9 13 < 10°. For 
simplicity we have further set 5 — 4>i — 4>2 = 0. 

Within the context of the Constrained Minimal 
Supersymmetric Standard Model (CMSSM), univer- 
sality of the soft SUSY breaking parameters is im- 
posed at a high-energy scale Mx, which we choose 
to be the SU(2) — U(l) gauge coupling unification 
scale (Mx ~ 2 x 10 16 GeV). Instead of scanning 
over the full CMSSM parameter space (generated 
by -M1/2, Mo, Aq, tan/3, sign/i), we considered spe- 
cific choices for the latter parameters, given by some 
of the "Snowmass Points and Slopes" (SPS) Q cases 
defined in Table [J 



Table I Values of M1/2, Mo, Ao, tan/3, and sign(/i) for the 
SPS points considered in the analysis. 



SPS 


M 1/2 (GeV) 


M (GeV) 


Ao (GeV) 


tan/3 




la 


250 


100 


-100 


10 


> 


lb 


400 


200 





30 


> 


2 


300 


1450 





10 


> 


3 


400 


90 





10 


> 


4 


300 


400 





50 


> 


5 


300 


150 


-1000 


5 


> 



Regarding our computation of the LFV observ- 
ables [lj, it is important to stress the following points: 

• It is a full one-loop computation of the branch- 
ing ratios (BRs), i.e., we include all contribut- 
ing one-loop diagrams with the SUSY particles 



flowing in the loops. For the case of L 



the analytical formulae can be found in [E llQ |. 
Regarding the lj —> 3k decays, the complete 
set of diagrams (including photon-penguin, Z- 
penguin, Higgs-penguin and box diagrams) and 
formulae are given in [6|. 

• The computation is performed in the physical 
basis for all SUSY particles entering in the loops. 
In other words, we do not use the Mass Insertion 
Approximation (MIA). 

• To obtain the low-energy parameters of the 
model the full renormalisation group equations 



(RGEs), including relevant terms and equations 
for the neutrinos and sneutrinos, are firstly run 
down from Mx to m^. At the seesaw scale 
(in particular at mN 3 ), we impose the bound- 
ary condition of Eq. (fTJ). After the decoupling 
of the heavy neutrinos and sneutrinos, the new 
RGEs are then run down from m,N 1 to the EW 
scale, at which the observables are computed. 
More concretely, we do not use the Leading Log 
Approximation (LLog), but rather numerically 
solve the full one- loop RGEs. 

• The numerical implementation of the above pro- 
cedure is achieved by means of the public For- 
tran code SPheno2.2.2 [ll|, which has been 
adapted in order to fully incorporate the right- 
handed neutrino (and sneutrino) sectors, as well 
as the full lepton flavour structure @. 

• The SPhcno code has been further enlarged by 
additional subroutines that compute the LFV 



branching ratios for all the lj — * kj and lj 



3k 



channels @. We have also included subrou- 
tines [l| to implement the requirement of suc- 
cessful baryogenesis (which we define as having 
ns/n 7 € [10 -10 , 10 -9 ]) via thermal leptogene- 
sis in the presence of upper bounds on the re- 
heat temperature, and to ensure compatibility 
with present bounds on lepton electric dipole 
moments (EDMs): EDM eur < (6.9x 1(T 28 , 3.7x 
10" 19 ,4.5 x 10" 17 )e.cm 

In what follows we present our main results for the 
case of hierarchical heavy neutrinos. We also include a 
comparison with present bounds on LFV rates 
El EE Ell and their future sensitivities [H,E1 
collected in Table ILTl 

Table II Present bounds and future sensitivities for the 
LFV processes. 



13, 
20, 



14, 
21] 



LFV process 


Present bound Future sensitivity 


BR(n -f e7) 
BR(r -> e-y) 
BR(r^/i 7 ) 
BR(/i -> 3 e) 
BR(r -> 3e) 
BR(r -> 3/x) 


1.2 x 10~ n 1.3 x 10" 13 
1.1 x 10" 7 10" 8 

6.8 x 10" 8 10" 8 
1.0 x 10~ 12 10" 13 
2.0 x 10" 7 10" 8 

1.9 x 10" 7 10" 8 



3. Results and Discussion 

Here we focus on the sensitivity of the BRs to #13, 
and on the dependence on other relevant parameters, 
which, for the case of hierarchical heavy neutrinos, 
are the heaviest mass r7ijv 3 , tan/3, 6\ and O2 (using 
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Figure 1: BR(/i — *• ef) and BR(/i — > 3e) as a function of 813 (in degrees), for SPS la (dots), lb (crosses), 2 (asterisks), 3 
(triangles), 4 (circles) and 5 (times). A dashed (dotted) horizontal line denotes the present experimental bound (future 
sensitivity) . 



the R parameterisation of [8|). The other input see- 
saw parameters mjvi, mN 2 an d $3: P ia Y a secondary 
role since the BRs do not strongly depend on them. 
Finally we comment on the hints on the SUSY seesaw 
parameters that can be derived from a measurement 
of the BRs and #13. 

For R = 1, the predictions of the BRs as func- 
tions of in the experimentally allowed range of #13, 
0° < 0i3 < 10° are illustrated in Fig. [TJ In this figure 
we also include the present and future experimental 
sensitivities for the channels. We clearly see that the 
BRs of fi — > e7 and /1 — > 3e are extremely sensitive to 
013, with their predicted rates varying many orders of 
magnitude along the explored 0i3 interval. The BRs 
of r — > e 7 and r — > 3 e channels are also sensitive 
to 0i3, but experimentally less challenging. The other 
LFV channels, r — > /17 and r — > 3/x, are nearly in- 
sensitive to this parameter (see [U). In the case of 
fi — > e7 this strong sensitivity was previously pointed 
out in Ref. [5(. In [6|, working within a full RGE 
approach, it was noticed that fi — > e7 and /1 — > 3e 
were the channels that, in addition to manifesting a 
clear 0i3 dependency, were the most promising from 
the experimental detection point of view. 

The most important conclusion from Fig. [1] is that, 
for this choice of parameters, the predicted BRs for 
both muon decay channels, fi —> ej and fi — > 3e, are 
clearly within the present experimental reach for sev- 
eral of the studied SPS points. The most stringent 
channel is manifestly \x — > ej where the predicted BRs 
for all the SPS points are clearly above the present 
experimental bound for 0i3 > 5°. With the expected 
improvement in the experimental sensitivity to this 
channel, this would happen for 0i3 > 1°. 

In addition to the small neutrino mass generation, 
the seesaw mechanism offers the interesting possibility 
of baryogenesis via leptogenesis [H • Thermal leptoge- 



nesis is an attractive and minimal mechanism to pro- 
duce a successful BAU, even compatible with present 
data, n B /n 7 w (6.10 ± 0.21) x 10~ 10 In the su- 
persymmetric version of the seesaw mechanism, it can 
be successfully implemented provided that the follow- 
ing conditions can be satisfied. Firstly, Big Bang Nu- 
cleosynthesis gravitino problems have to be avoided, 
which is possible, for instance, for sufficiently heavy 
gravitinos. Since we consider the gravitino mass as a 
free parameter, this condition can be easily achieved. 
In any case, further bounds on the reheat tempera- 
ture, Trh, still arise from decays of gravitinos into the 
lightest supersymmetric particle (LSP). In the case of 
heavy gravitinos and neutralino LSP masses in the 
range 100-150 GeV (which is the case of our work), 
one obtains Trh ^ 2 x 10 10 GeV. In the presence 
of these constraints on Trh, the favoured region by 
thermal leptogenesis corresponds to small (but non- 
vanishing) complex R- matrix angles 0^. For vanish- 
ing £/mns CP phases the constraints on R are basi- 
cally 1 02 1, 1 03 1 ;$ lrad (mod ir). Thermal leptogen- 
esis also constrains mjv x to be roughly in the range 
[10 9 GeV, 10 x T RH ] (see also [HEl). 

In [l| we have explicitly calculated the produced 
BAU in the presence of upper bounds on the re- 
heat temperature Tr,h- We have furthermore set as 
"favoured BAU values" those that are within the in- 
terval [10~ 10 , 10~ 9 ], which contains the WMAP value, 
and chosen the value of — 10 10 GeV in most of 
our analysis. Similar studies of the constraints from 
leptogenesis on LFV rates have been done in [251 ] . 

For very small values of m Vl (m Vl ~ 0(1O -5 eV)) a 
baryon asymmetry in the range 10~ 10 to 10~ 9 can be 
obtained for a considerable region of the |02| parame- 
ter space, with the BRs exhibiting a clear sensitivity 
to the value of 0i3 [l[ . On the other hand, the situa- 
tion changes dramatically for larger values of m Vl . 
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Figure 2: BR(/i — > e'y) as a function of \9i\, for arg #2 = 
{71-/8 , 7r/4 , 3-7r/8} (dots, times, diamonds, respectively) 
and #13 = 0°, 5° (blue/darker, green/lighter lines). We 
take m l , 1 — 10~ 3 eV. In all cases black dots represent 
points associated with a disfavoured BAU scenario and a 
dashed (dotted) horizontal line denotes the present exper- 
imental bound (future sensitivity). 



Figure 3: BR(/i — > ej) as a function of mjv 3 for SPS la, 
with m ui = 10 -5 eV and m vi = 10~ 3 eV (times, dots, 
respectively), and #13 = 0°, 5° (blue/darker, green/lighter 
lines). Baryogenesis is enabled by the choice 62 = 
0.05 e°- 2i (6»i = 3 = 0). On the upper horizontal axis we 
display the associated value of (1^)33. A dashed (dotted) 
horizontal line denotes the present experimental bound 
(future sensitivity). 



In Fig. we display the dependence of the most 
sensitive BR to f9 13 , BR(^ — > e-f), on \9 2 \. We con- 
sider two particular values of 6*13, 813 = 0° , 5° and 
choose SPS la. Motivated from the thermal leptogen- 
esis favoured # 2 -regions [J, we take < \0 2 \ < tt/4, 
with arg # 2 = {tt/8 , 7r/4 , 37r/8}. We choose = 
10 -3 eV, while for the heavy neutrino masses we take 
m N = (io 10 , 10 u , 10 14 ) GeV. 

While for smaller values of \0 2 \ the branching ratio 
displays a clear sensitivity to having #13 equal or dif- 
ferent from zero (a separation larger than two orders 
of magnitude for \0 2 \ < 0.05), the effect of r5 13 is di- 
luted for increasing values of |t9 2 |. For \0 2 \ ^ 0.3 the 
BR(/i — > 67) associated with 6*13 = 5° can be even 
smaller than for 6*13 = 0°. This implies that in this 
case, a potential measurement of BR(/i — * e 7) would 
not be sensitive to #13. Similar results were obtained 
for 63, but for shortness are not shown here. 

Concerning the EDMs, which are clearly non- 
vanishing in the presence of complex 0i, we have 
checked that all the predicted values for the electron, 
muon and tau EDMs are well below the experimental 
bounds. 

We now consider the dependence of BR(/i — » cy) 
on mN 3 - As displayed in Fig.[3J there is a strong sen- 
sitivity of the BRs to mjv 3 - In fact, the BRs vary by 
as much as six orders of magnitude in the explored 
range of 5 x 10 11 GeV < m Na < 5 x 10 14 GeV. Notice 
also that for the largest values of m/y 3 considered, the 
predicted rates for — > ej enter into the present ex- 
perimental reach. Although not shown here, it is also 
worth mentioning that by comparing our full results 
with the LLog predictions, we found that the LLog ap- 



proximation dramatically fails in some cases |l| . Sim- 
ilar effects were also noticed in (2(| [27j ■ 

Regarding the tan (3 dependence of the BRs we 
obtained that the BR grow as tan 2 (3. In fact, 
the hierarchy of the BR predictions for the sev- 
eral SPS points (as already manifest in Fig[T]) is 
dictated by the corresponding tan/3 value, with a 
secondary role being played by the given SUSY 
spectra. We found the following generic hierarchy: 
BRsps4 > BRspsib <L BRspsia > BRsps3 <L BRsps2 
> BRspss- 

Let us now address the question of whether a joint 
measurement of the BRs and #13 can shed some light 
on experimentally unreachable parameters, like mjv 3 . 
The expected improvement in the experimental sensi- 
tivity to the LFV ratios supports the possibility that 
a BR could be measured in the future, thus provid- 
ing the first experimental evidence for new physics, 
even before its discovery at the LHC. The prospects 
are especially encouraging regarding /1 -> e7, where 
the experimental sensitivity will improve by at least 
two orders of magnitude. Moreover, and given the 
impressive effort on experimental neutrino physics, a 
measurement of 613 will likely also occur in the fu- 
ture [H]. 

Given that, as previously emphasised, fx — > is 
very sensitive to 6*13 , whereas this is not the case for 
BR(r — > fij), and that both BRs display the same 
approximate behaviour with m/v 3 and tan/3, we have 
studied the correlation between these two observables. 
This optimises the impact of a 613 measurement, since 
it allows to minimise the uncertainty introduced from 
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not knowing tan /3 and mjv 3 , and at the same time of- 
fers a better illustration of the uncertainty associated 
with the i?-matrix angles. In this case, the correla- 
tion of the BRs with respect to ttin 3 means that, for a 
fixed set of parameters, varying mjv 3 implies that the 
predicted point (BR(t —tfij), BR(/i — > ej)) moves 
along a line with approximately constant slope in the 
BR(t — > /i7)-BR(/i — > ej) plane. On the other hand, 
varying #13 leads to a displacement of the point along 
the vertical axis. 

In Fig. [H we illustrate this correlation for SPS la, 
choosing distinct values of the heaviest neutrino mass, 
and scanning over the BAU-enabling i?-matrix angles 
(setting 63 to zero) as 

< |0i| < tt/4, -7J-/4 < arg6»! < tt/4, 
< |0 2 | < tt/4 , < arg0 2 < tt/4, 

m N3 = 10 12 , 10 13 , 10 14 GeV. (2) 

We considered the following values, #13 = 1°, 3°, 5° 
and 10°, and only included in the plot the BR predic- 
tions which allow for a favourable BAU. Other SPS 
points have also been considered but they are not 
shown here for brevity (see [l[). We clearly observe in 
Fig. 0] that for a fixed value of m n 3 , and for a given 
value of 6*i3, the dispersion arising from a #1 and 62 
variation produces a small area rather than a point in 
the BR(r — > /1 7)-BR(/j ej) plane. 

The dispersion along the BR(r — > /17) axis is of 
approximately one order of magnitude for all #13. In 
contrast, the dispersion along the BR(/i — > ej) axis 
increases with decreasing $13, ranging from an or- 
der of magnitude for #13 = 10°, to over three orders 
of magnitude for the case of small #13 (1°). From 
Fig. 2] we can also infer that other choices of mjv 3 
(for 0i3 £ [1°,10°]) would lead to BR predictions 
which would roughly lie within the diagonal lines de- 



picted in the plot. Comparing these predictions for 
the shaded areas along the expected diagonal "corri- 
dor" , with the allowed experimental region, allows to 
conclude about the impact of a 0i3 measurement on 
the allowed/excluded mjv 3 values. 

The most important conclusion from Fig. [4] is that 
for SPS la, and for the parameter space defined in 
Eq. @ , an hypothetical $13 measurement larger than 
1°, together with the present experimental bound on 
the BR(/i — > e-f), will have the impact of exclud- 
ing values of mN 3 si 10 14 GeV. Moreover, with the 
planned MEG sensitivity, the same $13 measurement 
can further constrain mjv 3 i$ 3 x 10 12 GeV. The im- 
pact of any other $13 measurement can be analogously 
extracted from Fig. [U 

As a final comment let us add that, remarkably, 
within a particular SUSY scenario and scanning over 
specific 0i and 02 BAU-enabling ranges for various 
values of 0i3, the comparison of the theoretical pre- 
dictions for BR(/i — > cy) and BR(r — > ^7) with 
the present experimental bounds allows to set 0i3- 
dependent upper bounds on mjv 3 - Together with the 
indirect lower bound arising from leptogenesis consid- 
erations, this clearly provides interesting hints on the 
value of the seesaw parameter TOjv 3 ■ With the planned 
future sensitivities, these bounds would further im- 
prove by approximately one order of magnitude. 

Ultimately, a joint measurement of the LFV branch- 
ing ratios, 0i3 and the sparticle spectrum would be a 
powerful tool for shedding some light on otherwise un- 
reachable SUSY seesaw parameters. It is clear from all 
this study that the interplay between LFV processes 
and future improvement in neutrino data is challeng- 
ing for the searches of new physics. 
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